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Abstract 


Strong  Planning  and  Forecast  Horizons 
in  a  Convex  Production  Planning  Problem 

by 

J.  T.  Teng,  G.  L.  Thompson,  and  S.  P.  Sethi 


\ 

We  consider  a  production-inventory  planning  problem  with  time 
varying  demands,  convex  production  costs,  and  a  warehouse  capacity 
constraint.  It  is  solved  by  use  of  the  Lagrangian  form  of  the 
maximum  principle.  The  possible  existence  of  strong  planning  and 
forecast  horizons  is  demonstrated.  When  they  exist,  they  permit 
the  breaking  up  of  the  whole  problem  into  a  set  of  smaller 
problems  which  can  be  solved  independently,  because  optimal 
decisions  up  to  a  strong  planning  horizon  are  completely 
independent  of  demand  data  beyond  the  next  forecast  horizon.  A 
forward  branch  and  bound  algorithm  is  developed  to  determine  such 
horizons  and  to  solve  the  whole  problem. 


Accession  For 

«TIS  GRA&I 
PTIC  TA15 
Un.-i!!v:'.;:icod 
Justification.  - 


la 


[> 


1 


0 

0 

0 

0 

0 

o 

o 

a 

| 

!  o 

; 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


1  .  Introduction 

This  paper  deals  with  a  production  planning  problem  which  has 
a  time-varying  demand.  In  particular  we  consider  the  problem  of 
finding  a  production  schedule  over  a  finite  period  [0,T]  for  a 
product  having:  (1)  deterministic  exogenous  demand,  (2)  strictly 
convex  production  costs,  (3)  strictly  increasing  inventory  costs 
and  (4)  and  warehouse  capacity  constraint.  This  is  an  extension 
of  a  problem  solved  by  Modigliani  and  Hohn  [53.  Kleindorfer  and 
Lieber  [1]  have  treated  a  similar  problem  using  the  extrapolation 
approach  in  an  optimal  control  framework;  see  also  Thompson  and 
Sethi  [11].  We  formulate  the  problem  as  an  optimal  control 
problem  and  solve  it  by  using  the  Lagrangian  form  of  the  maximum 
principle.  First,  we  characterize  the  form  of  the  optimal 
production  rate  rule  and  inventory  level  rule.  We  then  find 
strong  planning  and  strong  forecast  horizons  which  permit  the 
decomposition  of  the  whole  problem  into  a  set  of  smaller  problems, 
because  the  optimal  decisions  during  the  period  up  to  the  strong 
planning  horizon  are  completely  independent  of  the  data  beyond  the 
corresponding  forecast  horizon.  Finally,  we  derive  a  forward 
branch  and  bound  algorithm  to  solve  the  whole  problem  and  give  a 
numerical  example  to  illustrate  the  algorithm. 

Earlier  planning  horizon  results  for  production- inventor y 
problems,  obtained  by  using  dynamic  programming  and  variational 
arguments,  have  been  given  by  Kunreuther  and  Morton  [2,  33. 
Modigliani  and  Hohn  [53,  and  Wagner  and  Whitin  [123.  Other 
planning  horizon  results,  derived  by  applying  optimal  control 
theory,  were  given  by  Kleindorfer  and  Lieber  [13,  Lieber  [4], 
Morton  [6],  Pekelman  [7,  83,  and  Thompson,  Sethi,  and  Teng  [10]. 
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2.  The  Model  and  Its  Solution 

Consider  a  firm  producing  a  homogeneous  good  and  having  an 

inventory  warehouse.  Define  the  following  quantities: 

I(t):  inventory  level  at  time  t  (state  variable), 

u(t):  production  rate  at  time  t  (control  variable), 

d(t):  demand  rate  at  time  t  (exogeneous  variable);  bounded 
and  differentiable  for  all  t>0, 

T:  length  of  planning  period, 

f(u):  strictly  convex  nonnegative  increasing  production 

cost;  twice  differentiable  for  all  u>0, 

h(I):  strictly  increasing  nonnegative  inventory  holding 

cost;  differentiable  for  all  t>0, 

W:  upper  bound  of  its  warehouse  capacity;  W>0. 

Suppose  the  firm  wants  to  minimize  its  production  and  inventory 
costs  to  meet  the  given  exogeneous  demand  rate.  The  resulting 
problem  is  mathematically  equivalent  to  minimizing  the  following 
expression : 

Jsf  J[f(u(t))  +  h(I(t))]dt  (1) 

subject  to 

Ex ]  I ( t) =u( t) -d( t) ;  I(0)=I0<W,  (2) 

tPl]  I>0,  Cp 23  I<W,  [c]  u>0 ,  (3) 

where  a  dot  denotes  the  first  derivative  with  respect  to  time, 
x(t)  is  the  adjoint  variable  of  (2),  and  p  ^  <  t ) ,  p2(t)  and  c(t)  are 
the  Lagrange  variables  of  the  corresponding  constraints  as  given 
in  (3). 
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2.1  The  Necessary  Conditions  for  an  Optimal  Solution 

To  apply  the  maximum  principle,  see  C 9 3 ,  we  write  the 
Hamiltonian  function  of  (1)  as 

H=-f(  u) -h(  I  )+A  (u-d)  .  CM) 

The  Lagrangian  then  will  be 

L=H-p 1 (u-d)+p 2( d-u) +; u  .  (5) 

The  following  necessary  conditions  hold  (see  [9])  for  an  optimal 
solution : 


(a L/a  u)  =0  =  -f '  (  u)  +x  +p  ^ -p  2 ,  or 


g(A  +p 


if  x+p  1-p2<f'  (0) 
if  x  +p  1-P2>f*  (0) 


(6) 


where 


gsCfT1  , 


(7) 


which  it  may  be  noted  from  the  assumption  on  f,  is  a  strictly 
increasing  nonnegative  function;  the  adjoint  equation  satisfies, 


A  =-a  L/a I=h* (I) 


(8) 


and  the  transversality  conditions  are 


A (T) >0  and  A (T)I(T)=0; 


(9) 


the  complementarity  and  nonnegativity  conditions  are 

•  • 

p  1 1 1  p  ^  I  •  Pg^W-I),  p  2^  »  cu=0,  and  p  ^ c  >  0 ; 


r 

a  t 


(10) 


il 


/  the  adjoint  variable  x  is  continuous  except 
1  possibly  at  an  entry  or  an  exit  to  the  boundary 
<  conditions  l(t)=0  or  I(t)=W,  the  variables  I  and  (11) 

x+p.j-p2  are  continuous  everywhere,  and  moreover 
^  p  ^  <0  and  p  2^0 . 


Note  that  t.  is  an  initial  or  entry  time  to  I(t)=W  if 
1  w  — — * - *  — *" 

I(t.")<W  and  I(t.  )=I(t,+.)=W,  t^  is  a  final  or  exit  time  to 

1W  1W  1W  IW  -  - 

I(t)=W  if  I(t  ~)=I(t~  )  s'W  and  I(t-+)<W.  The  definitions  of  t.n 
fw  fw  fw  lU 

and  are  similar  and  are  obtained  by  changing  W  to  0. 

2.2  Optimal  Policies  for  Three  Possible  Cases 

There  are  only  three  different  possible  cases  for  the  values 
of  I(t)  for  some  time  interval  S.  We  shall  discuss  the  optimal 
policies  in  each  of  these  cases. 

Case  1  0<I(t)<W  for  all  tc  S 
This  implies  p1=p2  =  0  an^ 


r  u*  =  0  and  I  =  -d  if  x  <f' (0) 

lu*  =  g(x)  and  i  =  g(x)-d  if  x>f»(0) 


(12) 


Case  2  l(t)=0  and  l(t)=0  for  all  teS. 
This  implies 

p  2=0 ,  u*  =  d,  and  f'(d):X  +p,. 


(13) 
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J 
0 

Since  o  <0,  we  have  i+p^h'd).  Using  (13),  this  case  can  occur 
.only  if 

0 

£  fM (d) d<h ' ( I )  (14) 

0 

Case  3  I(t)=W  and  l(t)  =  0  for  all  te  S 
This  implies  that 


0 

J 
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0 

0 

0 
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p  ^  =0 ,  u  =d ,  and  f' (d)=i  -p  (15) 

From  the  fact  that  p'2<0,  we  have  the  following  conclusion: 

Case  3  can  happen  only  if 

f" ( d) d>h ' ( I )  . 

Simple  economic  interpretations  of  Cases  2  and  3  are  as 
follows:  In  Case  2,  f"(d)d<h’(I)  implies  f"(u)u<h'(0)  because 

u=d  and  1=0.  Taking  integral  of  both  sides  over  any  interval 
( t1 ,  tg) C  S  ,  we  have 

fc2  fc2 

f’  (u(t2))-f  (u(t1))=/t^f"(u)udt</tii  h’  (O)dtsh'  (OXtg-t^ 
or  1  1  (17) 

f'(u(t2))<f,(u(t1))+h'(0)(t2-t1) . 

In  (17),  f'(u(t2))  is  the  cost  of  making  an  additional  small 
quantity  at  time  t2  by  producing  it  at  time  t2-  On  the  other  hand 
f 1  ( u( t^ ) )  +  h' (0) ( t2«t ^ )  is  the  cost  of  making  an  additional  small 


0 


quantity  available  at  time  t2  by  producing  it  at  time  t^  and 
storing  it  until  time  t2,  hence  including  an  inventory  cost 
h'(0)(t2-t  ).  Therefore,  if  inequality  (17)  holds,  then  it  is  not 
worthwhile  to  build  up  inventory  for  the  future,  i.e.,  keeping  the 
warehouse  empty  is  the  optimal  strategy,  by  a  similar  argument, 
if  f" ( d) d>h' ( I )  holds  then  it  is  desirable  to  build  up  inventory. 
2.3  Theoretical  Results 

We  shall  characterize  the  optimal  trajectories  of  the 
problem.  First,  let  us  define 

IQs/ Jd(t)dt 

which  is  the  total  demand  during  the  planning  period. 

Theorem  1 .  The  following  three  statements  are  equivalent: 

(a)  i0>V 

(b)  I (T ) >0  and  x(T)<f’(0)  , 

( c)  u  ( t) =0  for  all  t . 

Proof .  We  prove  that  (a)  implies  (b)  only.  Other 
implications  are  trivial.  If  I(T)>0  then  by  (9)  we  have 
X (T)=0<f ' (0) .  If  I(T) =0  then  I0=Ig  and  u*=0  for  all  t.  Hence, 

X  (T)  <  f '  (0)  . 

Corollary  1  IQ<T0  implies  I(T)=0 
"Proof.  If  I(T)>0  then  x (T)=0< f' (0) .  From  Theorem  1,  we  have 
IQ> IQ ,  which  leads  to  a  contradiction . 

Theorem  2  We  can  conclude  that  I<0  for  all  te[0,  T]  under 


each  of  the  following  cases: 


(ii)  IQ<I  and  f"(d)d<h’(I)  for  all  te t 0 ,  T] 

Proof .  Case  (i)  immediately  follows  form  Theorem  1.  We  now 
prove  case  (ii).  Suppose  not,  i.e.,  I(t)>0  in  some  interval  (j, 
j+e)c[0,T].  This  can  happen  only  in  Case  1  of  Section  2.2. 

Hence,  I=g(x )-d>0  in  (j,  j+e  )  .  Since  f">0  from  the  strictly 
convexity  of  f  we  have  that  f'  is  a  strictly  increasing  function. 
Thus,  from  (7),  g(x )-d>0  if  and  only  if  x>f’(d).  Because 

x  =h'  (I)>d[  f'  (d)  ]/dt=f"(d)d  for  all  t  and  x>f'(d)  in  (j,j+e) 
we  get  x>f'(d)  for  all  t>j.  Therefore,  l(t)>0  for  all  t>j  which 
implies  I(T)>0  and  leads  to  a  contradiction  to  Corollary  1. 

Corollary  2  Under  the  cases  specified  in  Theorem  2,  the 
optimal  path  for  I(t)  is  character ized  by  one  of  the  following 
three  cases: 

(а)  I ( t) >0  for  all  t  (This  case  can  happen  if  and  only  if  Iq>Iq)  , 

(б)  I ( t )  =  0  for  all  t  (This  case  can  occur  if  and  only  if  1^=0), 

(c)  I ( t) >0  on  [ 0 , t  ^ )  and  l(t)=0  on  [t^.T]  (This  case  can  happen 

if  and  ony  if  0<Iq<Iq) 

Next,  we  shall  explore  the  case  in  which  f " ( d) d> h ’ ( I ) ,  i.e., 
the  case  in  which  it  pays  to  store  inventory.  For  convenience,  we 
define  the  function 

*  (t)s(x+p  1-p2)(t)  (18) 

Then  we  have  the  following  results. 

Theorem  3.  If  d  is  continuous  and  nondecreasing,  then  we 

have : 

(1)  *  is  continuous  and  nondecreasing, 
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(2)  u  is  continuous  and  nondecreasing. 

Proof.  If  0<I<W  for  all  te  S  then  by  (12)  we  know  that  iji  =x  is 
an  increasing  function  because  x=h'(I)>0.  Otherwise,  1=0  or  I=W 
implies  that  =  f *  ( d )  is  also  a  nondecreasing  function  by  (13)  and 
(15).  Since  i|<  must  be  continuous  everywhere,  see  (11),  we  have 
that  *  is  continuous  and  nondecreasing. 

By  (6)  and  (7),  we  know  that  a  continuous  and  nondecreasing  ^ 

* 

will  imply  a  continuous  and  nondecreasing  u  because  g  is  a 
strictly  function. 

Corollary  3  If  d  is  continuous  and  nondecreasing,  then  the 

» 

optimal  control  path  for  u  ( t)  satisfies  one  of  the  following 
three  cases: 

(I)  u  =0  for  all  t 

(II)  u*> 0  for  all  t 

(III)  u  =0  for  t<tQ  and  u  >0  for  t>t^ 

Proof .  It  immdeiately  follows  from  Theorem  3- 

Theorem  4  If  f"(d)d>h'(I)  for  all  te[0,T]  then  we  have  the 

following  results:  (1)  l(t)>0  almost  everywhere  on  [0,T],  (2) 

« 

u  (t)  will  be  one  of  the  three  cases  as  shown  in  Corollary  3- 
Proof.  f"(d)d>h'(I)  implies  Case  2  in  Section  2.2  cannot 

happen.  Therefore,  l(t)>0  almost  everywhere  on  [0,T].  Again, 

•  •  # 

f"(d)d>h'(I)  implies  d>0.  Hence  u  (t)  satisfies  one  of  the  three 
cases  shown  in  Corollary  3- 

From  Corollary  2  and  Theorem  4,  we  can  easily  get  the 

following  results: 

Theorem  5  Assume  Iq<Iq  an^  t’  is  choosen  so  that 
f"(d)d>h’(I)  for  t<t’  and  f"( d) d<h' (I)  for  t>t' ,  i.e.,  it  pays  to 
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'  9 

store  production  before  but  not  after  t' .  Then  the  optimal  path 
for  I(t)  satisfies  one  of  the  following  two  cases: 

(A)  l(t)>0  for  all  te[0,T] 

( B)  There  exists  t^>t’  such  that  l(t)>0  on  and  l(t)=0  on 

C 1 2 ,  T  3  . 

Theorem  6  Assume  Iq<Iq  and  t'  is  choosen  so  that 
f"(d)d<h’  (I)  for  t<t'  and  f"(d)d>h»(I)  for  t>t' ,  i.e.  it  pays  to 
store  production  after  but  not  before  t' .  Then  the  optimal  path 
for  I ( t )  satisfies  one  of  the  following  three  cases: 

(a)  l(t)>0  for  all  teCO.T], 

(B)  There  exists  t^<t  ’  such  that  I ( t) =0  on  [ 0  , t^ 3  and  I ( t) >0  on 
(t3,  T), 

(c)  There  exist  t4<t^<t'  such  that  l(t)>0  on  [0,t^),  l(t)=0  on 
[tjj,  t_],  and  l(t)>0  on  (t^,T). 

In  general,  the  interval  [0,T]  will  contain  many  subintervals 

•  • 

in  which  f"( d) d<h' (I)  or  f"  (  d)  d>  h' ( I )  .  By  repeatedly  applying  the 
results  of  Theorems  5  and  6  we  can  construct  the  solution  by 
piecing  together  different  solutions  obtained  from  applications  of 
these  theorems. 
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3.  Strong  Planning  Horizon  Theorem 

» 

If  optimal  production  decisions  during  [0,t  ]  can  be  shown  to 

««  « 

be  completely  independent  of  the  data  beyond  t  >t  ,  then  we  call 

•  *  * 

t  a  strong  planning  horizon  and  t  a  strong  forecast  horizon. 


In  this  section,  we  will  obtain  the  strong  planning  and 
strong  forecast  horizons  for  the  problem. 

t  *  * 

Lemma  1.  Let  t  and  t  be  the  two  consecutive  entries  at 
which  inventory  hits  the  boundary  constraint.  Furthermore,  let 

ft  JHt  ft  ft  ft  ft  ft  ft 

I(t  )=0  and  I(t  )=W,  or  I(t  )=W  and  I ( t  )=0.  Suppose  t  <t 

ft  ft  ft 

Then  t  is  a  strong  planning  horizon,  and  t  is  a  strong  forecast 


horizon . 

*  *# 

Proof .  We  will  prove  the  case  in  which  I(t  )=W  and  I(t  )=0. 
By  using  an  analogous  argument  we  could  prove  the  other  case  in 

ft  ft  ft 

which  I(t  )=0  and  I(t  )=W.  We  may  assume,  without  loss  of 

* 

generality,  that  t  is  the  first  time  at  which  the  process  enters, 

i.e.,  satisfies,  the  boundary  condition  I(t)=W.  Let  the  optimal 

value  of  x (0)  for  the  original  problem  be  denoted  by  xQ.  If  we 

t  ,  * 

can  prove  that  x ( t) =x Q+; Qh ' ( I ) dt  on  (0,t  )  is  still  an  optimal 

solution  to  any  problem  having  the  same  demand  rate  information  as 

,  #» 

the  original  problem  on  [0,t  ]  and  regardless  of  its  value  after 

** 

t  ,  then  we  are  done.  If  the  new  optimal  value  of  x  (0)  were  xQ 

~  ~  * 

with  Xq<Xq,  then  the  new  inventory  level  I<I  on  [0,t  ]  because 

^  ft  ^  ft 

>4<  =  > u  >u  by  using  (6)  and  (7) 

•  *  '»  i 

=  >I  =  u  -d> u  -d=I . 

ft 

Hence,  I(t  )<W.  For  an  illustration  see  Figure  1.  Since 

» 

* =x -p 2=  f 1  ( d)  on  (t  ,tQ)  where  tQ  is  the  time  of  exit  from  the 
inventory  constraint  I(t)=W,  we  know 
♦  =  f"(d)d>h*  (I)>h>  (I)=*  . 


Figure  1 


•  i  -  » 

This  implies  that  * >* ,  I>I,  and  I>I  on  (t  ,t0).  Similarly,  we  can 
get 


I < I ,  I< I ,  and  I<W  on  [  t, 


ft* 

,t  ] 


This  implies 

>  ,»  -  t**r  t«*  . 

I(t  )-I(t  )=/  Idt<;  Idt=I(t  )-I(tn)=-w 

U  V  “  t  U 

«  •» 

or  I < t  )<0  and  leads  to  a  contradiction.  By  using  an  analogous 
argument  we  can  prove  that  if  the  new  optimal  value  of  x (0)  were 


xQ  with  x >x  then  it  leads  to  a  contradiction  again. 

Theorem  7  (Strong  Planning  Horizon  Theorem) 

Let  t,  and  t-  be  the  initial  and  final  times  to  I(t)=W 
iw  fw 

respectively,  i.e.  I(t)=W  for  te [ t .  ,t«  ].  Suppose  t.-  is  the 
next  time  at  which  the  inventory  reaches  l(t)=0  boundary.  Then 
all  t  before  tfw  are  strong  planning  horizons  and  all  t  after  t^ 
are  strong  forecast  horizons,  i.e.  t^w  is  a  maximal  strong 
planning  horizon  and  tiQ  is  a  minimal  forecast  horizon.  The 
Theorem  also  holds  when  0  and  W  are  interchanged. 

Proof.  If  t  „  is  not  a  strong  planning  horizon  then  there 

I  w 

exists  a  new  optimal  solution  such  that  toE^tiw,tfw^  new 

exit  form  I(t)=W  boundary.  The  reasons  are  as  follows:  to<tiw 
will  contradict  to  Lemma  1  and  tQ>tfw  implies  that  t^w  is  a  strong 
planning  horizon.  Similarly,  using  the  same  arguments  as  in  Lemma 
1  we  can  prove  that  the  new  exit  tQ  will  lead  to  a  contradiction. 
4.  A  Forward  Branch  and  Bound  Algorithm 

The  solution  is  trivial  if  Iq>Iq.  Therefore,  we  may  assume, 
without  loss  of  generality,  that  Iq<Iq.  Suppose  j  is  a  feasible 
solution  to  problem  (1)  during  [0,t^].  Then  it  is  easy  to  show 
(proof  omitted)  that  a  lower  bound  for  this  policy  j  during  [0,T] 
is 


i(J) 


Cf(u(  j)) 


h(I(j))]dt  + 


f(u)[T-t.], 


(19) 


where 


u=[/ l  ddt  -  I(ti)]/(T-ti) 


f 


(20) 


13 


u(j}  and  are  the  production  rate  and  inventory  level 

corresponding  to  the  policy  j,  respectively.  From  (19)  and  (20), 

we  also  have  that  J0=f(uQ)T  is  a  lower  bound  to  the  value  of  J 

—  • 

where  Uq= ( Iq-Iq) /T .  On  the  other  hand,  u  =d  for  all  t  is  a 

T 

feasible  solution  to  problem  (1).  Thus,  Jn=h(I-) tn+/  f(d)dt 

U  0  U  t-Q 

^0 

where  tn  is  the  solution  of  / -  d  dt=In,  is  an  upper  bound  of  J. 
We  are  now  in  a  position  to  present  theualgorithm.  For 

convenience,  we  shall  say  that  the  vertex  v.  is  fathomed  if  and 

J  - 

only  if  no  further  exploration  from  this  vertex  can  be  profitable 

Otherwise,  we  shall  say  that  v.  is  unfathomed  or  alive. 

J  -  - 


Forward  Branch  and  Bound  Algorithm 


Step  0  (Initialization)  Begin  at  the  live  vertex  vQ,  where 
J=JQ  and  J  =  Jq.  Go  to  Step  1. 

Step  1  (Branching)  Assume  that  t^  is  the  first  entry  time  to  the 
constraint  I(t)=W  or  l(t)=0. 

Case  1.1  Suppose  that  t^  is  the  entry  time  of  I=W.  Then  t^ 
should  satisfy  the  following  equations  and  contraints: 


x(t"i)=f»(d(ti)) 


*  (t)=x  (tT)-/t.1h’ (I)dt  for  all  t<t. 

1  w  —  1 


W-I0>/J[g(x)-d]dt>-I0  for  all  t<tt 


L 


3 

3 

9 

3 

9 

3 

9 

9 

9 

9 

9 

o 

3 

9 

9 

9 

9 

9 

9 

9 

9 

9 
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W-I„=/  i[g(x)-d]dt  (24) 

u  0 

f"(d(t  ))d>h'(W)  at  t-tt  (25) 

If  there  exists  one  or  more  solutions,  then  we  keep  each  such 
solution  as  successor  vertices  of  v^  and  go  to  Case  1.2. 

Otherwise,  there  are  no  solutions  for  this  branch,  terminate  this 
branch  and  go  to  Case  1.2. 

Case  1  .2  Suppose  that  t^  is  the  entry  of  1  =  0.  Then  t^ 
satisfies  the  same  constraints  as  in  Case  1.1,  except  that  (24) 
and  (25)  must  be  replaced  by  (26)  and  (27),  respectively. 


I0=/0i(gU)-d]dt 

(26) 

f" ( d) d<h' (0)  at  t=t.*T 

(27) 

Again,  we  keep  all  solutions  as  successor  vertices  of  v^ ,  if  any, 
and  go  to  Step  3. 

Step  2  (Fathoming  by  Bound)  Check  every  new  live  vertex  v  ^  .  If 
as  in  (19)  is  larger  than  or  equal  to  J,  then  the  vertex  v. 

J 

l-3  fathomed ,  i.e.,  no  further  exploration  from  this  vertex  can  be 
profitable.  Go  to  Step  4. 

Step  3  (Update  Bound)  Check  each  new  live  vertex  v..  If  t^^sT 

- c -  J  i 

and  then  the  vertex  v^  is  fathomed.  If  t^^=T  and  J^^<J 

then  let  J=J^.  Go  to  Step  2. 


Step  4  (Branching)  If  no  live  vertices  exist,  go  to  Step  7; 
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) 

1 


> 

otherwise  select  a  live  vertex 

which  I  reaches  to  W  or  0.  If 

v^.  Let  t[^  be  the 
I(t^)=W  then  go  to 

last  time 

Step  5; 

at 

J 

otherwise  go  to  Step  6. 

3 

Step  5  Assume  that  t^w  is  the 

exit  time  from 

I(t)=W 

,  and  t^ 

is 

the  next  entry  time  of  1=0  or 

I=W.  As  in  Step 

1 ,  we 

have  two 

3 

cases . 

3 

Case  5.1  Suppose  t^  is  the  entry  time  to 

1  =  0, 

and  t_  <t 
fw 

i  ‘ 

3 

Solving  the  following  constraints,  we  can  find 

and  t . . 

l 

the  values  of 

lfw 

J 

) 

x(tfw+)=f.(d(tfw)) 

(28) 

) 

x(0  =  f’(d(t.)) 

(29) 

3 

x (t)-x (tfw+)=/^  h'(I)dt 
fw 

fcr  all  te ( tfw , 

t.  ) 
l 

(30) 

V 

3 

-W</^  [  g(x  )-d]  dt<0 

<w 

for  all  te (tfw, 

V 

(3D 

3 

f"(d)d>h'  (W) 

for  all  tc  (tpJ 

’  tfw) 

(32) 

3 

3 

3 

-Ws/ .  1 C g(x  )-d]dt 

(33) 

r 

3 

f”(d)d<h' (0) 

at  t=t^lT 

(31*) 

*  If  above  constraints  have  one  or  more  solutions,  then  we  save  them 

Q 


I  o 

0 


f) 


3 

0 

0 

3 

9 


as  successor  vertices  of  v^  and  go  to  Case  5.2.  Otherwise, 
terminate  this  branch  and  go  to  Case  5.2. 

Case  5 .2  Let  t^  be  the  entry  time  to  I(t)=W.  This  case  is 
similar  to  Case  5.1,  except  that  (33)  and  (34)  are  replaced  by 
(35)  and  (36),  respectively. 


0=/ 


[g(x )-d]dt 


(35) 


o 

3 

3 

0 

0 

0 

o 

o 

o 

o 

(5 

3 


f"(d)d>h’  (W) 


at  t=t. 

l 


(36) 


If  there  exist  some  solutions  to  (tfw,  tp  in  Case  5.1  or  Case 

5.2,  then  we  keep  them  as  successor  vertices  of  v.  and  to  to  Step 

J 

3.  Otherwise,  v.  is  fathomed  and  go  to  Step  4. 

J 


Step  6  Assume  that  t _  is  the  exit  time  from  I(t)=0.  and  t.  is 

- f  q  x 

the  next  entry  time  to  1=0  or  I=W. 

Case  6.1  Suppose  that  t^  is  the  entry  time  to  1=0.  Solving 
the  following  constraints,  we  may  get  the  values  of  t  and  t^ . 

x(tf0+)=f'(d(tf0))  (37) 


x(t')=f’(d(t.)) 


(38) 


x  (t)-x  (tfQ+)=/ 


h' < I )dt 


for  all  tc  ( tfQ,  t^ 


(39) 


O 
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0 </*  [g(X  )-d]dt<W 
X0 

f"(d)d<h' (0) 

Or/.1  [ g(X  )-d]dt 

tfrt 


f"(d)d<h*  (0)  at  trt^T  (43) 

We  keep  all  solutions  as  successors  of  v  ,  if  any,  and  go  to  Case 

J 

6.2. 

Case  6 .2  Suppose  that  ^  is  the  entry  time  to  I(t)=W. 

Again,  this  case  is  similar  to  Case  6.1,  except  that  we  replace 
(42)  and  (43)  by  (44)  and  (45),  respectively. 


If  there  exist  any  solutions  to  (t^Q,  t^)  in  Case  6.1  or  Case  6.2 
then  save  them  and  go  to  Step  3*  Otherwise,  check  whether  (41)  is 
satisfied  by  setting  tf0  =  T  or  not.  If  yes,  let  I(t)r0  on 
(t^\  T)  be  a  feasible  solution  and  go  to  Step  3-  If  not,  v^  is 
fathomed  and  go  to  Step  4. 

Step  7  (Termination)  JrJ  which  is  optimal. 
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To  illustrate  that  the  algorithm  can  find  an  optimal 
solution,  we  will  solve  a  simple  numerical  example. 

An  Example 

Suppose  that  d< t) =60+10t-t2 ,  T=10,  f(u)=u2/2,  h(I)=I, 
W=125/6,  and  IQ  =  1 7 -  We  then  have  that  f'(u)  =  u,  f"(u)=1,  g(u)  =  u, 

X  =  1 ,  d= 1 0-2t ,  f"(d)d>h’  (I)  for  t<4.5,  and  f"  ( d)  d<  h*  ( I )  for  t>4.5. 
Step  0  J»f(15)* 10=1,125  and  JQ=29,027 

Step  1  Solving  Case  1.1,  we  have 

*(0  =  f'(d<t.)  )  =60-*-l  0t  . 

Substituting  the  result  into  (22)  and  (24),  we  obtain 


W-I0=/o1[(9trti)2-9t+t2]dt, 


and  t^=1.  Let  the  vertex  corresponding  this  case  is  v1 .  Next,  we 
find  that  there  is  not  feasible  solution  to  Case  1.2.  and  go  to 


Step  3. 


Step  3  t)  =1$T  and  go  to  Step  2. 

Step  2  J^<J  and  go  to  Step  4. 

Step  4  There  is  only  one  live  vertex  v^  ,  and  I(tP  )=W.  Go 


to  Step  5. 


Solving  (28)-(34)  simultaneously,  yields  t^w=2  and 


t^=7.  Let  v^  be  the  vertex  corresponding  to  th: s  case.  Again,  we 
find  that  there  are  no  feasible  solutions  to  Case  5.2.  Go  to  Step 


Step  3  t^=74T  and  go  to  Step  2. 

( 

Step  2  J  <J  and  go  to  Step  4  again. 

(2) 

Step  4  There  exists  only  one  live  vertex  v2,  and  I(t^  )=0. 


o 

Q 

O 

© 

o 

o 

o 

o 


o 

o 

© 

o 

o 

o 

o 

o 

o 

© 

o 

o 
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Go  to  Step  6. 

Step  6  There  are  no  feasible  solutions  to  both  Case  6.1  and 
6.2.  Since  f"(d)d<h’(0)  for  all  tc[7,  10].  We  know  that  l(t)=0 
on  [7,10]  is  a  feasible  solution,  let  be  the  corresponding 
vertex  to  this,  and  go  to  Step  3. 


Step  3  J=/i0[f(u*)+h(I*)  ]dt, 


where 


and 


i*  = 


0 

68+t 

74+t 

60+1 0t-t‘ 


if  t<1 
if  2<t<7 
otherwise , 


(46) 


(47) 


» 

I  =  < 


f  17+8t-9t2/2+t3/3 
125/6 

49/6+l4t-9t2/2+t3/3 

to 


if  t<1 
if  1 <  t<2 
if  2<t<7 
if  7< t< 1 0 . 


(48) 


Go  to  Step  2. 

Step  2  =  J  so  that  is  fathomed.  Go  to  Step  4. 

Step  4  There  is  no  live  vertex.  Go  to  Step  7. 

Step  7  J  as  in  (46)  is  the  optimal  solution.  This  shows 
that  t=2  is  a  maximal  strong  planning  horizon  and  t=7  is  a  minimal 
strong  forecast  horizon,  see  Figure  2. 
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5.  Conclusions 

In  this  paper  we  have  studied  a  general  production  inventory 
model  with  convex  costs  and  an  inventory  upper  bound.  We 
characterized  the  optimal  trajectories  and  showed  that  there  could 
be  strong  planning  and  forecast  horizons.  When  such  horizons 
exist,  we  showed  that  the  problem  can  be  decomposed  into  a  set  of 
smaller  problems.  In  Section  5  we  presented  a  forward  branch  and 
bound  algorithm  which  carries  out  this  decomposition.  The 
algorithm  is  illustrated  with  a  simple  example. 

It  would  be  possible  to  extend  the  results  of  this  paper  in 
several  different  ways.  For  instance  more  general  demand 
functions  could  be  considered.  Also,  if  the  state  equation  (2)  is 
replaced  by  I=u-d-kl,  where  k>0  is  the  rate  of  spoilage  of  the 
inventory,  our  strong  planning  horizon  theorem  is  still  true.  The 
other  theorems  in  this  paper,  with  suitable  modification,  are  also 
true . 


0 


0 

0 
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